Unraveling the Thermodynamics of the Folding and Interconversion of Human Telomere G-Quadruplexes.

Why human telomere DNA fragments fold into different G-quadruplex structures with parallel, hybrid, and antiparallel strand orientations depending on the temperature and concentration of co-solutes remains poorly understood. Similarly, the formation of intermediate structures along the folding or interconversion pathways is not well understood. Herein, we address these questions by introducing a conceptual framework, based on the global thermodynamic analysis of DSC and CD spectroscopy data, which led to a detailed description of the topological phase space (phase diagram) of the stability of the human telomere fragment 5'-AGGG(TTAGGG)3 -3' (Tel22). This framework clarifies the driving forces of quadruplex folding and interconversion processes over a wide range of temperatures and ion (K(+) , Na(+) ) and polyethylene glycol (PEG) concentrations and demonstrates their linkage to the human telomere DNA structural features.

Thermodynamic fingerprints of ligand binding to human telomeric G-quadruplexes.

Thermodynamic studies of ligand binding to human telomere (ht) DNA quadruplexes, as a rule, neglect the involvement of various ht-DNA conformations in the binding process. Therefore, the thermodynamic driving forces and the mechanisms of ht-DNA G-quadruplex-ligand recognition remain poorly understood. In this work we characterize thermodynamically and structurally binding of netropsin (Net), dibenzotetraaza[14]annulene derivatives (DP77, DP78), cationic porphyrin (TMPyP4) and two bisquinolinium ligands (Phen-DC3, 360A-Br) to the ht-DNA fragment (Tel22) AGGG(TTAGGG)3 using isothermal titration calorimetry, CD and fluorescence spectroscopy, gel electrophoresis and molecular modeling. By global thermodynamic analysis of experimental data we show that the driving forces characterized by contributions of specific interactions, changes in solvation and conformation differ significantly for binding of ligands with low quadruplex selectivity over duplexes (Net, DP77, DP78, TMPyP4; KTel22 ≈ KdsDNA). These contributions are in accordance with the observed structural features (changes) and suggest that upon binding Net, DP77, DP78 and TMPyP4 select hybrid-1 and/or hybrid-2 conformation while Phen-DC3 and 360A-Br induce the transition of hybrid-1 and hybrid-2 to the structure with characteristics of antiparallel or hybrid-3 type conformation.

Dominant Driving Forces in Human Telomere Quadruplex Binding-Induced Structural Alterations.

Recently various pathways of human telomere (ht) DNA folding into G-quadruplexes and of ligand binding to these structures have been proposed. However, the key issue as to the nature of forces driving the folding and recognition processes remains unanswered. In this study, structural changes of 22-mer ht-DNA fragment (Tel22), induced by binding of ions (K(+), Na(+)) and specific bisquinolinium ligands, were monitored by calorimetric and spectroscopic methods and by gel electrophoresis. Using the global model analysis of a wide variety of experimental data, we were able to characterize the thermodynamic forces that govern the formation of stable Tel22 G-quadruplexes, folding intermediates, and ligand-quadruplex complexes, and then predict Tel22 behavior in aqueous solutions as a function of temperature, salt concentration, and ligand concentration. On the basis of the above, we believe that our work sets the framework for better understanding the heterogeneity of ht-DNA folding and binding pathways, and its structural polymorphism.

Energetic basis of uncoupling folding from binding for an intrinsically disordered protein.

Intrinsically disordered proteins (IDPs) are proteins that lack a unique three-dimensional structure in their native state. Many have, however, been found to fold into a defined structure when interacting with specific binding partners. The energetic implications of such behavior have been widely discussed, yet experimental thermodynamic data is scarce. We present here a thorough thermodynamic and structural study of the binding of an IDP (antitoxin CcdA) to its molecular target (gyrase poison CcdB). We show that the binding-coupled folding of CcdA is driven by a combination of specific intramolecular interactions that favor the final folded structure and a less specific set of intermolecular contacts that provide a desolvation entropy boost. The folded structure of the bound IDP appears to be defined largely by its own amino acid sequence, with the binding partner functioning more as a facilitator than a mold to conform to. On the other hand, specific intermolecular interactions do increase the binding affinity up to the picomolar range. Overall, this study shows how an IDP can achieve very strong and structurally well-defined binding and it provides significant insight into the molecular forces that enable such binding properties.

Energetic basis of human telomeric DNA folding into G-quadruplex structures.

Recent theoretical studies performed on the folding/unfolding mechanism of the model telomeric human DNA, 5'-AGGGTTAGGGTTAGGGTTAGGG-3' (Tel22), have indicated that in the presence of K(+) ions Tel22 folds into two hybrid G-quadruplex structures characterized by one double and two reversal TTA loops arranged in a different way. They predicted a new unfolding pathway from the initial mixture of hybrid G-quadruplexes via the corresponding intermediate triplex structures into the final, fully unfolded state. Significantly, no experimental evidence supporting the suggested pathway has been reported. In the current work, we performed a comprehensive global thermodynamic analysis of calorimetric (DSC, ITC) and spectroscopic (CD) data obtained on monitoring the folding/unfolding of Tel22 induced by changes of temperature and K(+) concentration. We show that unfolding of Tel22 may be described as a monomolecular equilibrium three-state process that involves thermodynamically distinguishable folded (F), intermediate (I), and unfolded (U) state. Considering that calorimetric methods cannot distinguish between energetically similar G-quadruplex or triplex conformations predicted by the theoretical model one can conclude that our results represent the first experimental support of the suggested unfolding/folding mechanism of Tel22. This conclusion is confirmed by the fact that the estimated number of K(+) ions released upon each unfolding step in our thermodynamic model agrees well with the corresponding values predicted by the theoretical model and that the observed changes in enthalpy, entropy, and heat capacity accompanying the F → I and I → U transitions can be reasonably explained only if the intermediate state I is considered to be a triplex structural conformation.

Recognition of human tumor necrosis factor α (TNF-α) by therapeutic antibody fragment: energetics and structural features.

Human tumor necrosis factor α (TNF-α) exists in its functional state as a homotrimeric protein and is involved in inflammation processes and immune response of a human organism. Overproduction of TNF-α results in the development of chronic autoimmune diseases that can be successfully treated by inhibitors such as monoclonal antibodies. However, the nature of antibody-TNF-α recognition remains elusive due to insufficient understanding of its molecular driving forces. Therefore, we studied the energetics of binding of a therapeutic antibody fragment (Fab) to the native and non-native forms of TNF-α by employing calorimetric and spectroscopic methods. Global thermodynamic analysis of data obtained from the corresponding binding and urea-induced denaturation experiments has been supported by structural modeling. We demonstrate that the observed high affinity binding of Fab to TNF-α is an enthalpy-driven process due mainly to specific noncovalent interactions taking place at the TNF-α-Fab binding interface. It is coupled to entropically unfavorable conformational changes and accompanied by entropically favorable solvation contributions. Moreover, the three-state model analysis of TNF-α unfolding shows that at physiological concentrations, TNF-α may exist not only as a biologically active trimer but also as an inactive monomer. It further suggests that even small changes of TNF-α concentration could have a considerable effect on the TNF-α activity. We believe that this study sets the energetic basis for understanding of TNF-α inhibition by antibodies and its unfolding linked with the concentration-dependent activity regulation.

Influence of a hairpin loop on the thermodynamic stability of a DNA oligomer.

DSC was used to evaluate the mechanism of the thermally induced unfolding of the single-stranded hairpin HP = 5'-CGGAATTCCGTCTCCGGAATTCCG-3' and its core duplex D (5'-CGGAATTCCG-3')(2). The DSC melting experiments performed at several salt concentrations were successfully described for HP and D in terms of a three-state transition model HP↔I (intermediate state) ↔ S (unfolded single-stranded state) and two state transition model D↔2S, respectively. Comparison of the model-based thermodynamic parameters obtained for each HP and D transition shows that in unfolding of HP only the HP↔I transition is affected by the TCTC loop. This observation suggests that in the intermediate state its TCTC loop part exhibits significantly more flexible structure than in the folded state while its duplex part remains pretty much unchanged.

Kinetically governed polymorphism of d(G4T4G3) quadruplexes in K+ solutions.

It has been generally recognized that understanding the molecular basis of some important cellular processes is hampered by the lack of knowledge of forces that drive spontaneous formation/disruption of G-quadruplex structures in guanine-rich DNA sequences. According to numerous biophysical and structural studies G-quadruplexes may occur in the presence of K(+) and Na(+) ions as polymorphic structures formed in kinetically governed processes. The reported kinetic models suggested to describe this polymorphism should be considered inappropriate since, as a rule, they include bimolecular single-step associations characterized by negative activation energies. In contrast, our approach in studying polymorphic behavior of G-quadruplexes is based on model mechanisms that involve only elementary folding/unfolding transitions and structural conversion steps that are characterized by positive activation energies. Here, we are investigating a complex polymorphism of d(G(4)T(4)G(3)) quadruplexes in K(+) solutions. On the basis of DSC, circular dichroism and UV spectroscopy and polyacrylamide gel electrophoresis experiments we propose a kinetic model that successfully describes the observed thermally induced conformational transitions of d(G(4)T(4)G(3)) quadruplexes in terms of single-step reactions that involve besides single strands also one tetramolecular and three bimolecular quadruplex structures.

Model-based thermodynamic analysis of reversible unfolding processes.

Folding and unfolding of many biological macromolecules can be characterized thermodynamically, yielding a wealth of information about the stability of various conformations and the interactions that hold them together. The relevant thermodynamic parameters are usually obtained by employing spectroscopic and/or calorimetric techniques and fitting an appropriate thermodynamic model to the experimental data. In this work, we compare the traditional approach of fitting the thermodynamic model to experimental data obtained from each experiment individually and the global approach of simultaneously fitting the model to all available data from different experiments. On the basis of several specific examples of DNA and protein unfolding, we demonstrate that piece-by-piece verification of the proposed thermodynamic model using individual fits is frequently inappropriate and can result in an incorrect mechanism and thermodynamics of the studied unfolding process. We find that while the two approaches are complementary in some aspects of analysis global fitting is essential for the appropriate selection and critical evaluation of the model mechanism. Only a good global fit thus gives us confidence that the obtained thermodynamic parameters of unfolding have real physical meaning.

Energetics of MazG unfolding in correlation with its structural features.

MazG is a homodimeric alpha-helical protein that belongs to the superfamily of all-alpha NTP pyrophosphatases. Its function has been connected to the regulation of the toxin-antitoxin module mazEF, implicated in programmed growth arrest/cell death of Escherichia coli cells under conditions of amino acid starvation. The goal of the first detailed biophysical study of a member of the all-alpha NTP pyrophosphatase superfamily, presented here, is to improve molecular understanding of the unfolding of this type of proteins. Thermal unfolding of MazG monitored by differential scanning calorimetry, circular dichroism spectroscopy, and fluorimetry at neutral pH in the presence of a reducing agent (dithiothreitol) can be successfully described as a reversible four-state transition between a dimeric native state, two dimeric intermediate states, and a monomeric denatured state. The first intermediate state appears to have a structure similar to that of the native state while the final thermally denatured monomeric state is not fully unfolded and contains a significant fraction of residual alpha-helical structure. In the absence of dithiothreitol, disulfide cross-linking causes misfolding of MazG that appears to be responsible for the formation of multimeric aggregates. MazG is most stable at pH 7-8, while at pH <6, it exists in a molten-globule-like state. The thermodynamic parameters characterizing each step of MazG denaturation transition obtained by global fitting of the four-state model to differential scanning calorimetry, circular dichroism, and fluorimetry temperature profiles are in agreement with the observed structural characteristics of the MazG conformational states and their assumed functional role.

Driving forces of gyrase recognition by the addiction toxin CcdB.

Gyrase, an essential bacterial topoisomerase, is the target of several antibiotics (e.g. quinolones) as well as of bacterial toxin CcdB. This toxin, encoded by Escherichia coli toxin-antitoxin module ccd, poisons gyrase by causing inhibition of both transcription and replication. Because the molecular driving forces of gyrase unfolding and CcdB-gyrase binding were unknown, the nature of the CcdB-gyrase recognition remained elusive. Therefore, we performed a detailed thermodynamic analysis of CcdB binding to several fragments of gyrase A subunit (GyrA) that contain the CcdB-binding site. Binding of CcdB to the shorter fragments was studied directly by isothermal titration calorimetry. Its binding to the longer GyrA59 fragment in solution is kinetically limited and was therefore investigated via urea induced unfolding of the GyrA59-CcdB complex and unbound GyrA59 and CcdB, monitored by circular dichroism spectroscopy. Model analysis of experimental data, in combination with the relevant structural information, indicates that CcdB binding to gyrase is an enthalpic process driven mainly by specific interactions between CcdB and the highly stable dimerization domain of the GyrA. The dissection of binding energetics indicates that CcdB-gyrase recognition is accompanied by opening of the tower and catalytic domain of GyrA. Such extensive structural rearrangements appear to be crucial driving forces for the functioning of the ccd toxin-antitoxin module.

Diverse polymorphism of G-quadruplexes as a kinetic phenomenon.

Knowledge of forces that drive conformational transitions of G-quadruplexes is crucial for understanding the molecular basis of several key cellular processes. It can only be acquired by combining structural, thermodynamic and kinetic information. Existing biophysical and structural evidences on polymorphism of intermolecular G-quadruplexes have shown that the formation of a number of these structures is a kinetically controlled process. Reported kinetic models that have been used to describe the association of single strands into quadruplex structures seem to be inappropriate since the corresponding model-predicted activation energies turn out to be negative. By contrast, we propose here a novel kinetic model that successfully describes experimentally monitored folding/unfolding transitions of G-quadruplexes and gives positive activation energies for all elementary steps, including those describing association of two single strands into bimolecular quadruplex structures. It is based on a combined thermodynamic and kinetic investigation of polymorphic behavior of bimolecular G-quadruplexes formed from d(G4T4G4) and d(G4T4G3) strands in the presence of Na(+) ions, monitored by spectroscopic (UV, CD) and calorimetric (DSC) techniques. According to our experiment and model analysis the topology of the measured G-quadruplexes is clearly flexible with the conformational forms that respond to the rate of temperature change at which global unfolding/folding transitions occur.

What drives the binding of minor groove-directed ligands to DNA hairpins?

Understanding the molecular basis of ligand-DNA-binding events, and its application to the rational design of novel drugs, requires knowledge of the structural features and forces that drive the corresponding recognition processes. Existing structural evidence on DNA complexation with classical minor groove-directed ligands and the corresponding studies of binding energetics have suggested that this type of binding can be described as a rigid-body association. In contrast, we show here that the binding-coupled conformational changes may be crucial for the interpretation of DNA (hairpin) association with a classical minor groove binder (netropsin). We found that, although the hairpin form is the only accessible state of ligand-free DNA, its association with the ligand may lead to its transition into a duplex conformation. It appears that formation of the fully ligated duplex from the ligand-free hairpin, occurring via two pathways, is enthalpically driven and accompanied by a significant contribution of the hydrophobic effect. Our thermodynamic and structure-based analysis, together with corresponding theoretical studies, shows that none of the predicted binding steps can be considered as a rigid-body association. In this light we anticipate our thermodynamic approach to be the basis of more sophisticated nucleic acid recognition mechanisms, which take into account the dynamic nature of both the nucleic acid and the ligand molecule.

Energetics in correlation with structural features: the case of micellization.

Understanding micellization processes at the molecular level has direct relevance for biological self-assembly, folding, and association processes. As such, it requires complete characterization of the micellization thermodynamics, including its correlation with the corresponding structural features. In this context, micellization of a series of model non-ionic surfactants (poly(ethylene glycol) monooctyl ethers, C(8)E(gamma)) was studied by isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC). The corresponding structural properties of C(8)E(gamma) micelles were investigated by small-angle X-ray scattering (SAXS). The C(8)E(gamma) micellization, characterized independently from ITC, DSC, and structural data, reveals that deltaH(M)(o) > 0, deltaS(M)(o) > 0, and deltaC(P)(M)(o) < 0, while the dissection of its energetics shows that it is primarily governed by the transfer of 20-30 C(8) alkyl chains from aqueous solution into the nonpolar core (r approximately 1.3 nm) of the spherical micelle. Moreover, thermodynamic parameters of micellization, estimated from the structural features related to the changes in solvent-accessible surface areas upon micellization, are in a good agreement with the corresponding parameters obtained from the analysis of ITC and DSC data. We have shown that the contributions to deltaS(M)(o) other than from hydration (deltaS(M)(other)(o)), estimated from experimental data, appear to be small (deltaS(M)(other)(o) < 0.1 deltaS(M)(other)(o)) and agree well with the theoretical estimates expressed as a sum of the corresponding translational, conformational, and size contributions. These deltaS(M)(other)(o) contributions are much less unfavorable than those estimated for a rigid-body association, which indicates the dynamic nature of the C(8)E(gamma) micellar aggregates. the dynamic nature of the C8EY micellar aggregates.

Erythropoietin unfolding: thermodynamics and its correlation with structural features.

Human erythropoietin (EPO) is a glycoprotein hormone considered to be the principal regulator of red blood cell formation. Although its recombinant version (rEPO) has been widely used for treatment of various anemias and its biological effects are relatively well-known, we know little about its biophysical properties and their relation to its structure. To gain a fuller understanding of the structural and functional properties of rEPO on the molecular level we followed its thermal and urea-induced unfolding at different pH (3.1-9.4) and urea concentrations (0-8 M) using spectropolarimetry, UV absorption, intrinsic emission fluorescence, and differential scanning calorimetry. Our results show that under a variety of conditions rEPO undergoes thermal or urea-induced denaturation that may be considered as a reversible two-state process characterized by unusually high (thermal) or moderate (urea-induced) extent of the residual structure. The highest thermal stability of the protein observed in aqueous solutions at physiological pH appears to be due to the largest difference in the extent of structure in the denatured and native state at this pH. The comparison between experimentally determined energetics of rEPO denaturation and its structure-based calculations indicates that the parametrization of thermodynamic quantities in terms of changes in solvent accessible nonpolar and polar surface areas resulting from protein unfolding can be successfully used provided that these changes are estimated from combination of experimentally determined deltaC(o)p and deltaH(o) values and not calculated from the structure of the protein's folded and assumingly fully unfolded state.

Energetics of structural transitions of the addiction antitoxin MazE: is a programmed bacterial cell death dependent on the intrinsically flexible nature of the antitoxins?

The Escherichia coli mazEF addiction module plays a crucial role in the cell death program that is triggered under various stress conditions. It codes for the toxin MazF and the antitoxin MazE, which interferes with the lethal action of the toxin. To better understand the role of various conformations of MazE in bacterial life, its order-disorder transitions were monitored by differential scanning calorimetry, spectropolarimetry, and fluorimetry. The changes in spectral and thermodynamic properties accompanying MazE dimer denaturation can be described in terms of a compensating reversible process of the partial folding of the unstructured C-terminal half (high mean net charge, low mean hydrophobicity) and monomerization coupled with the partial unfolding of the structured N-terminal half (low mean net charge, high mean hydrophobicity). At pH

Energetic diversity of DNA minor-groove recognition by small molecules displayed through some model ligand-DNA systems.

Energetics of interactions occurring in the model ligand-DNA systems constituted from distamycin A (DST), netropsin (NET) and the oligomeric duplexes d(GCAAGTTGCGATATACG)d(CGTATATCGCAACTTGC)=D#1 and d(GCAAGTTGCGAAAAACG)d(CGTTTTTCGCAACTTGC)=D#2 was studied by spectropolarimetry, UV-absorption spectroscopy and isothermal titration calorimetry. Model analysis of the measured signals was applied to describe individual and competitive binding in terms of populations of various species in the solution. Our results reveal several unprecedented ligand-DNA binding features. DST binds to the neighboring 5'-AAGTT-3' and 5'-ATATA-3' sites of D#1 statistically in a 2:1 binding mode. By contrast, its association to D#2 appears to be a 2:1 binding event only at the DST/D#2 molar ratios between 0 and 2 while its further binding to D#2 may be considered as a step-by-step binding to the unoccupied 5'-AAAAA-3' sites resulting first in DST3D#2 and finally in DST4D#2 complex formation. Competition between DST and NET binding shows that for the most part DST displaces NET from its complexes with D#1 and D#2. In contrast to the obligatory 1:1 binding of DST to the ligand-free 5'-AAAAA-3' sites observed at DST/5'-AAAAA-3' <1 the displacement of NET bound to the 5'-AAAAA-3' sites by added DST occurs even at the smallest additions of DST in a 2:1 manner. NET can also displace DST molecules but only those bound monomerically to the 5'-AAAAA-3' sites of DST3D#2. Actually, only half of these molecules can be displaced due to the simultaneous rebinding of the displaced DST to the unreacted 5'-AAAAA-3' sites in DST3D#2. Binding of DST and NET to D#1 and D#2 is an enthalpy driven process accompanied by large unfavorable (DST), small (NET) or large favorable (NET binding to 5'-AAAAA-3') entropy contributions and negative deltaCP degrees that are reasonably close to deltaCP degrees predicted from the calculated changes in solvent-accessible surface areas that accompany complex formation. Although various modes of DST and NET binding within D#1 and D#2 are characterized by significant energetic differences they seem to be governed by the same driving forces; the hydrophobic transfer of ligand from the solution into the duplex binding site and the accompanying specific non-covalent ligand-DNA and ligand-ligand interactions occurring within the DNA minor groove.